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ABSTRACT 



Neutrino oscillations affect light element synthesis through the ^-process in 
supernova explosions. The 7 Li and 11 B yields produced in a supernova explo- 
sion of a 16.2 M & star model increase by factors of 1.9 and 1.3 in the case of 
large mixing angle solution with normal mass hierarchy and sin 2 2^13 >2x 10~ 3 
compared with those without the oscillations. In the case of inverted mass hi- 
erarchy or nonadiabatic 13-mixing resonance, the increment of their yields is 
much smaller. Neutrino oscillations raise the reaction rates of charged- current 
^-process reactions in the region outside oxygen-rich layers. The number ratio 
of 7 Li/ n B could be a tracer of normal mass hierarchy and relatively large #13, 
still satisfying sin 2 26*13 < 0.1, through future precise observations in stars having 
strong supernova component. 
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Subject headings: neutrinos — nuclear reactions, nucleosynthesis, abundances - 
supernovae: general 

1. Introduction 

In the final evolutionary stage of massive stars, most region of the stars except col- 
lapsing core explodes as supernova explosions. The collapsing core releases its gravitational 
energy with gigantic amount of neutrinos. The emitted neutrinos interact with nuclei in the 
exploding material and new species of nuclei are produced; this synthetic process is called 
the //-process (Domogatsky, Eramzhyan, & Nadyozhin 1978; Woosley et al. 1990). There are 
several species produced through the //-process. For light elements, 7 Li and 11 B are mainly 
produced through the //-process (Woosley et al. 1990; Woosley & Weaver 1995; Yoshida, 
Emori, & Nakazawa 2000; Rauscher et al. 2002; Yoshida et al. 2004; Yoshida, Kajino, & 
Hartmann 2005). Some 19 F is also produced through the //-process (Woosley et al. 1990; 
Woosley & Weaver 1995; Rauscher et al. 2002). For neutron-deficient heavy nuclei 138 La and 
180 Ta are also produced through charged- current interactions with v e (Goriely et al. 2001; 
Heger et al. 2005). Neutrino-driven winds from proto- neutron stars are considered to be one 
of the promoting sites for r-process heavy elements (e.g., Woosley et al. 1994; Takahashi, 
Witti, & Janka 1994; Otsuki et al. 2000; Terasawa et al. 2004). 

Supernova explosion is one of the important sites for supplying 7 Li and n B as well as 
Galactic cosmic rays, AGB stars, and novae during Galactic chemical evolution (GCE) (e.g., 
Fields et al. 2000). In previous studies, we showed that the amounts of 7 Li and n B strongly 
depend on the neutrino energy spectra and the total neutrino energy (Yoshida et al. 2004, 
2005). We also constrained the neutrino energy spectra from the gravitational energy of 
a proto-neutron star and GCE models (Yoshida et al. 2005). In these studies it has been 
assumed that the neutrino spectra do not change in the supernova ejecta. 

On the other hand, recent remarkable progress in neutrino experiments has confirmed 
the phenomenon of neutrino oscillations (e.g., McKeown & Vogel 2004). The experiments 
on atmospheric neutrinos (e.g., Ashie et al. 2004), solar neutrinos (e.g., Ahmed et al. 2004), 
and reactor neutrinos (e.g, Apollonio et al. 2003; Araki et al. 2005) constrained most of 
parameter values in neutrino oscillations, such as squared mass differences and the mixing 
angles. Resultantly, the large mixing angle (LMA) solution turns out to be a unique solution 
for the 12- and 23-mixings. However, mass hierarchy between 1 and 3 mass eigenstates 
has not been clarified (e.g., Ashie et al. 2004) and only upper limit of sin 2 2^13 has been 
determined (Apollonio et al. 2003). 
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Supernova neutrino is another promoting target for neutrino experiments. When SN 
1987A occurred, Kamiokande group and 1MB group found eleven and eight events of neutrino 
detection (Hirata et al. 1987; Bionta et al. 1987). Owing to the development of neutrino 
experiments, much larger events of the neutrino detection are expected when a supernova 
explosion occurs in neighboring galaxies. In order to evaluate the neutrino flux and their 
energy dependence, neutrino oscillations in supernova explosions have been investigated 
qualitatively (Dighe & Smirnov 2000) and quantitatively (Takahashi et al. 2001). They 
showed that in the case of adiabatic resonance for sin 2 29 13 (LMA-L in Takahashi et al. 
2001) the transition probability from v e to v^ T changes from to almost 1 in the O/C layer 
in their supernova model. Finally, the energy spectrum of u e changes to the one close to 
the v^ T spectrum emitted from neutrino sphere. In the case of nonadiabatic resonance for 
sin 2 2#i3 (LMA-S in Takahashi et al. 2001), the change of neutrino spectra is smaller. The 
effects from mass hierarchy and from the change of the density profile due to the shock 
propagation were also investigated (Takahashi & Sato 2003; Takahashi et al. 2003). 

Since neutrino oscillations change neutrino spectra, it is expected that the amounts of 7 Li 
and n B change by the effect of neutrino oscillations (Yoshida et al. 2006). During supernova 
explosions, neutral-current reactions such as 4 He(z/, u'p) 3 B. and 4 He(z/, z/'n) 3 He are important 
for 7 Li and n B production (Yoshida et al. 2004). We note here that the total reaction 
rates of neutral-current reactions do not change by the neutrino oscillations. The energy 
spectrum summed up in all neutrinos and antineutrinos does not change by the oscillations. 
On the other hand, the reaction rates of charged-current reactions such as 4 He(z/ e , e~p) i Hc 
and 4 He(z/ e , e + n) 3 H are expected to increase by neutrino oscillations. As shown in Dighe & 
Smirnov (2000) and Takahashi et al. (2001), the mean energies of v e and v e increase by the 
neutrino oscillations. This increase will raise the efficiency of the 7 Li and n B production. 
If we obtain some clear signals of neutrino oscillations in the abundances of 7 Li and n B, 
we would constrain the parameter values of neutrino oscillations from observations of light 
elements. This is a new procedure to constrain neutrino oscillation parameters completely 
different from the detections of supernova neutrinos. 

In the present study, we investigate light element synthesis in supernova explosions 
taking account of the change of the neutrino spectra due to neutrino oscillations. We also 
evaluate the dependence of the yields of 7 Li and n B in the supernova ejecta on the mixing 
angle sin 2 2#i 3 and mass hierarchy. 

We set the luminosity and the energy spectrum of each flavor of neutrinos just emitted 
from a proto-neutron star in §2. We also set the parameter values of neutrino oscillations 
from the results of recent neutrino experiments. We explain a supernova explosion model 
and a nuclear reaction network for light element synthesis. We mention the cross sections of 
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charged- current reactions of the //-process to evaluate the reaction rates including neutrino 
oscillations. In §3, we show the transition probabilities of neutrino flavors with different 
values of sin 2 29 13 and mass hierarchy. We also discuss the effect of the oscillations on the 
reaction rates of charged-current //-process reactions. In §4, we show the calculated mass 
fraction distribution of 7 Li and n B taking account of neutrino oscillations. Then, we show 
the dependence of the 7 Li and n B yields on sin 2 29 13 and mass hierarchy. We also show the 
dependence on the temperatures of v e and P e just emitted from proto-neutron star. In §5, 
we discuss the neutrino oscillations with supernova shock propagation and show the change 
of the 7 Li and n B yields by this effect. We also discuss the 7 Li and n B yields related to 
observations of stars which have traces of supernova explosions and supernova remnants. 
Finally, we conclude our study in §6. 



2. Supernova Model and Parameters 

2.1. Models of Supernova Neutrinos Emitted from a Proto-Neutron Star 

We use a model of supernova neutrinos just emitted from a proto-neutron star based 
on the models in previous studies on the //-process nucleosynthesis (Yoshida et al. 2004, 
2005, 2006). We here set up the energy spectra of three flavors of neutrinos when they just 
emitted from the proto-neutron star before they have been affected by neutrino oscillations in 
passing through the envelope. The neutrino luminosity is assumed to decrease exponentially 
with the decay time of r v = 3 s (after Woosley et al. 1990). The total neutrino energy E v 
is set to be 3.0 x 10 53 ergs which is almost equal to the gravitational binding energy of a 
1.4 M® neutron star (e.g., Lattimer & Prakash 2001). The neutrino luminosity is equally 
partitioned for each flavor of neutrinos. The neutrino energy spectra are assumed to obey 
Fermi-Dirac (FD) distributions with zero-chemical potentials, similarly to previous studies 
on the //-process. The influence of nonzero chemical potentials has been discussed in Yoshida 
et al. (2005). The temperature of z/ MjT and z/ MjT , T v T , is set to be 6.0 MeV. This neutrino 
temperature is adopted so that the production of n B from supernovae satisfies appropriate 
for GCE constraint on the light elements (Yoshida et al. 2005). The temperatures of v e 
and z/ e , T Ue and T Pe , are set to be 3.2 MeV and 5.0 MeV, with which we have investigated 
light element synthesis and the r-process heavy element synthesis (Yoshida et al. 2004). We 
investigate detailed influences of neutrino oscillations using this set of neutrino temperatures. 

We also use the temperature of 4.0 MeV for v e and z/ e , which has been adopted in 
Woosley & Weaver (1995), Rauscher et al. (2002), and others. It has been indicated that 
the cooling of proto-neutron stars makes the temperatures of each flavor of neutrinos closer 
to each other (e.g., Keil et al. 2003). Therefore, we consider the case of T„ e and T Pe to be 4.0 
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MeV and 5.0 MeV, respectively, too. For comparison, we set T Ve and T Pe as 3.2 MeV and 
4.0 MeV, corresponding to the case where only T Ve is changed. 



2.2. Parameters for Neutrino Oscillations 



In the present study, we consider the change of the neutrino spectra due to 3-flavor 
neutrino oscillations 



the 



dx 









( 








\ 


V , 


H 







Am^c 4 







V T ) 






V° 





Am§ lC 4 
2e v 


I 



I ±V2G F (h C y 



V 



m u 



\ 

j 



u = 



C12C13 S12C13 

— s 12 c 23 ~~ Ci2S23 s 13 c 12 c 23 ~~ s 12 s 23 s 13 
S12S23 — c 12 c 23 s 13 ~ C12S23 — •Sl2C23 s 13 




(1) 

(2) 



where Am-, = rri} 



m?, rrii is the mass of i-eigenstate neutrinos, e u is the neutrino en- 
ergy, Gp is Fermi constant, p is the density, Y e is electron fraction, m u is the atomic mass 
unit, Sij = sin0y, and = cos^-. Positive and negative signs in the potential term 
±a/2G f(^c) 3 pY e / 'm u correspond to the cases of neutrinos and antineutrinos. In this for- 
mulation the squared mass differences and the mixing angles are parameters. The values 
of these parameters except #13 and the sign of Am^ have been precisely determined from 
recent neutrino experiments by Super-Kamiokande (Ashie et al. 2004), SNO (Ahmed et al. 
2004), and KamLAND (Araki et al. 2005). We assume the mass differences and the mixing 
angles as follows: 



Am 2 21 



7.9 x 10" 5 eV 2 and 



Am 2 31 



±2.4 x 10" 



eV^ 



(3) 



and 



sin 2 2#i2 = 0.816 and sin 2 2# 23 = 1.0. 



(4) 



This parameter set corresponds to the LMA solution of neutrino oscillations. The sign of 
Am 2 x has not been determined from recent neutrino experiments. Thus, we consider both 
normal mass hierarchy and inverted mass hierarchy. The positive value and the negative 
value of Am^ correspond to normal mass hierarchy, i.e., mi < m 2 < m 3 , and inverted mass 
hierarchy, i.e., < mi < m,2, respectively. From CHOOZ experiment (Apollonio et al. 
2003) only upper limit of sin 2 26 ) 13 has been determined. We investigate the influence of 
changing 13 in the range of 

< sin 2 2#i3 < 0.1. (5) 
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2.3. Supernova Explosion Model 

We calculate light element nucleosynthesis using the same supernova explosion model 
adopted in our previous studies (Yoshida et al. 2004, 2005). The presupernova model is 14E1 
model, which is a 16.2 M & star just before the supernova explosion and corresponds to SN 
1987A (Shigeyama & Nomoto 1990). The shock propagation of the supernova is calculated 
using a piecewise parabolic method code (Colella & Woodward 1984; Shigeyama et al. 1992). 
The explosion energy is set to be 1 x 10 51 ergs. The mass cut is located at 1.61 M . Note 
that it is not needed to calculate the structure inside the mass cut. In the present study, we 
calculate neutrino oscillations with the density structure of the presupernova model. When 
the shock front is in inner high density region, there is no influence of the shock propagation 
to the neutrino oscillations. This is because the oscillation amplitude is too small despite 
the presence of the shock wave in such a high density region. Its effect will be discussed in 
detail in §5. 



2.4. Nuclear Reaction Network for the //-Process with Neutrino Oscillations 

Our nuclear reaction network consists of 291 species of nuclei from n, p, to Ge. The 
included nuclear species and their associated nuclear reactions are the same as those in 
Yoshida et al. (2004), except charged-current //-process reactions on 4 He and 12 C, which will 
be discussed later in this subsection. The rates of neutral-current //-process reactions and 
those of the other charged-current reactions are adopted from the rates with the assumption 
of FD distribution of the neutrino spectra. As mentioned in Introduction, the rates of the 
neutral-current reactions do not change by neutrino oscillations. 

When we consider neutrino oscillations, the energy spectra of v e and V e emitted from 
a proto-neutron star change from FD distributions and charged-current //-process reaction 
rates are not merely simple functions of neutrino temperatures. In such a condition, we need 
to evaluate a neutrino flux as a function of the neutrino energy. The number flux of i-flavor 
neutrinos (i = e, /i, r) emitted from the proto-neutron star with the energy of e v can be 
written as 

d4n = L Vi 1 el 

de u 47rr 2 F 3 { Vu J{kT u ^ exp(^ - Vut ) + 1 1 ' 

and 

F ^ = L exp(/-t) + l' <7) 
where L v . is the luminosity of i-flavor neutrinos, r is the radius, k is the Boltzmann constant, 
T Vi is the neutrino temperature, and r\ Vi is the degeneracy factor (e.g., Balantekin & Yiiksel 
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2005). Here we define Pij(r; e v ) as the oscillation probability from i-flavor to j-flavor at the 
radius of r and the energy of e v . Taking the cross section of a //-process reaction as a function 
of £„, cr Ue (e u ), we evaluate the reaction rate with neutrino oscillations X Ue as 

^r 2 FMikT^ Jo exp(^ - Vui ) + 1 U 

In the present study, we consider the effect on neutrino oscillations for charged-current reac- 
tions of 4 He, 4 He(z/ e , e~p) 3 He and 4 He(z/ e , e + n) 3 H, and the reactions of 12 C, 12 C(z/ e , e~p) u C, 
12 C(z/ e ,e- 7 ) 12 N, 12 C(z/ e ,e + n) n B, and 12 C(u e , e+ 7 ) 12 B. However, the detailed data of the 
cross sections of these reactions as functions of e v have not been reported and one can find 
only energy averaged reaction rates as functions of neutrino temperature. Therefore, we 
adopt analytical approximation to the cross sections with respect to e u as 

/ x f a(e v - e th ) b for e v > e th , , , 

°"^ = {o for s„<e, h , (9) 

where e t h is the threshold energy of each reaction. Coefficients are determined so as to fit the 
rate of the corresponding reaction with the assumption of the FD distributions to the values 
tabulated in the 1992 work by R. D. Hoffman & S. E. Woosley, HW92 K The coefficients 
and the thresholds are listed in Table 1. When we evaluate the reaction rates from the cross 
sections with assumptions of the FD distributions, these values are in reasonable agreement 
with those of HW92 within ±9 %. 



3. Neutrino Oscillations in Presupernovae 

In order to evaluate the reaction rates of the charged-current //-process reactions with 
neutrino oscillations, we numerically solve the oscillations of the neutrinos passing through 
a presupernova star in Runge-Kutta method and based on analytical expression of Kimura, 
Takamura, & Yokomakura (2002a,b). We use the density profile of 14E1 model of Shigeyama 
& Nomoto (1990). The density profile is shown in Fig. 1 in Shigeyama & Nomoto (1990). 

When neutrinos pass through the stellar interior, there are resonances of neutrino os- 
cillations where the transition probabilities change largely. The resonance density p res is 
determined from the squared mass difference and the neutrino energy. The resonance den- 
sity is written as 

m u Am 2 i c 4 cos26» i7 - „ / Am|\ / 1 MeV\ , , . 

p ms Y e = ^ = 6.55xl0 6 3 4)i cos20 i7 - gem" 3 . 10 



1 http://www-phys.llnl.gov/Research/RRSN/nu_csbr/neu_rate.html. 
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Electron number density n e relates to the density and electron fraction through n e = pY e /m u . 
The transition probabilities depend on the adiabaticity of the resonance strongly. The adi- 
abaticity is estimated using the adiabaticity parameter 7. The flip probability Pf, which 
means the probability that a neutrino in one mass eigenstate changes to another mass eigen- 
state is written as (e.g., Dighe & Smirnov 2000) 

P f = exp 

and 

_ Am^c 3 sin 2 2% 
2he u cos 29ij 

When 7 ^> 1, i.e., the flip probability Pf is very small (Pf <C 1), the resonance is adiabatic. 
When the flip probability is close to unity, the resonance is nonadiabatic. 

Figure 1 shows the transition probability of v e to z/ e , P ee , the sum of the transition 
probabilities of — > v e and v T — > z/ e , P xc , and the corresponding transition probabilities for 
antineutrinos P ee b and P xe b in the case of e v = 50 MeV and normal mass hierarchy. There 
are resonances of Am^ and #13, which is called H-resonance, at p res = 6.3 x 10 2 g cm~ 3 
and of Am 2 x and 6> 12 , which is called L-resonance, at p rcs = 8.9 g cm" 3 . The H-resonance 
is located in innermost region of the He/C layer (~ 3.8M ) and the L-resonance is in the 
He/C layer (~ 5.4M ). Both of H- and L-resonances appear for neutrinos and do not 
appear for antineutrinos. The adiabaticity of H-resonance depends on the value of sin 2 29 13 . 
L-resonance is always adiabatic in our parameter set. We see in the mass coordinate range 
of M r < 3.8M that the amplitude of neutrino oscillation is very small and that the flavor 
exchange practically does not occur for both neutrinos and antineutrinos. 

In the mass coordinate region of M r > 3.8M , the characteristics of the transition 
probabilities depend on mass hierarchy and the adiabaticity of H-resonance. In the case of 
sin 2 2^13 = 1 x 10~ 2 (see Fig. la), H-resonance appears for neutrinos and is adiabatic. The 
transition from v e to v^ T occurs completely, i.e., P ec becomes almost zero in the He/C layer. 
At the same time, the transition probability from or v T to v e becomes large, i.e., P xe is close 
to 1. On the other hand, there is no resonance for antineutrinos. The transition probabilities 
P eeb and P xc b gradually change in the He/C layer, and about 30% of antineutrinos change 
flavors. 

In the case of sin 2 2#i 3 = 1 x 10~ 6 (see Fig. 16), H-resonance is nonadiabatic. So, 
complete change from u e to T does not occur. The transition probabilities for neutrinos 
change gradually as a function of the mass coordinate similarly to antineutrinos. Finally, 
about 70% of neutrinos change flavors. The transition probabilities for antineutrinos are the 
same as the case of sin 2 2#i 3 = 1 x 10 -2 . 



(11) 



/ l/n e )(dn e /dr) 



(12) 
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Figure 2 shows the transition probabilities in the case of inverted mass hierarchy. In 
this case, H- resonance appears for antineutrinos. When sin 2 2# 13 is equal to 1 x 10~ 2 (see 
Fig. 2a), H-resonance is adiabatic and almost all v e change to and v T . The transition 
probability P xc b also becomes close to unity. On the other hand, neutrinos change their 
flavors gradually in the He-layer because of no appearance of H-resonance. When sin 2 1Q\z 
is equal to 1 x 10~ 6 , H-resonance is nonadiabatic. The change of the transition probabilities 
with increasing the mass coordinate is the same as that in the case of normal mass hierarchy 
and the nonadiabatic H-resonance. In the case of nonadiabatic resonance, the flavor change 
occurs as if there is no resonance. 

The increase in charged-current //-process reaction rates strongly depends on the adi- 
abaticity of H-resonance. If H-resonance is adiabatic, the reaction rates of charged-current 
reactions for u e (normal mass hierarchy) and v e (inverted mass hierarchy) become much 
larger than those without the oscillations in the He-layer. If H-resonance is nonadiabatic, 
the increase in the rates of corresponding charged-current reactions would be much smaller. 
Therefore, the adiabaticity of H-resonance affects the final yields of 7 Li and n B. 

We note that the shock propagation effect on neutrino oscillations in supernova ejecta 
is not considered in this study. We do not consider that the change of neutrino oscillations 
due to the shock propagation would change the yields of 7 Li and n B. When neutrinos pass 
through the O/Ne layer, the amplitude of neutrino oscillations is very small because the 
density profile, where the density is much larger than H-resonance density, does not affect 
neutrino oscillations. Thus, when the shock wave moves inside the O/C layer, the shock 
wave does not affect the transition probabilities of neutrinos and antineutrinos even if we 
consider the shock propagation. After the shock wave arrives at the O/C layer, the shock 
wave will change the transition probabilities. We set the decay time of the neutrino flux 
to be 3 s which is to be compared with the shock arrival time to the O/C layer about 5 
s. Because of this time lag, more than 80% of neutrinos pass through the supernova ejecta 
before the shock arrives at the O/C layer. We will discuss the influence of the shock wave 
on neutrino oscillations in the supernova in §5. 



4. Yields of 7 Li and n B 

4.1. Mass Fraction Distributions of 7 Li and n B 

We show the change of the mass fractions of 7 Li and n B due to neutrino oscillations. 
Figure 3 shows the mass fraction distributions of 7 Li and n B. Panels (a) and (b) correspond 
to the case of normal mass hierarchy, and panels (c) and (d) correspond to inverted mass 
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hierarchy. In these figures we show the mass fraction distributions of 7 Li and its isobar 'Be 
separately. The mass fraction distributions of n B and n C are also drawn separately. 

Let us first discuss the case of normal mass hierarchy and sin 2 2#i 3 = 1 x 10~ 2 (adiabatic 
H-resonance for neutrinos) shown by thick lines in Figs. 3a and 36. We find the increase 
in the mass fractions of 7 Li, 7 Be, n B, and n C compared with those without the neutrino 
oscillations in the He layer. The mass fraction of 7 Li with the neutrino oscillations is larger 
by about a factor of 1.2 than that without the oscillations in most region of the He layer (see 
Fig. 3a). In the range of 4.3M Q < M r < 4.6M Q , the increment degree is larger. In the O-rich 
layers of M r < 3.8M Q , we do not see any clear differences due to neutrino oscillations in the 
mass fraction. The transition probabilities to other neutrino flavors are very small because 
of high density in this region (see Fig. 1). The obtained yield of 7 Li is 1.79 x 1O~ 7 M . 
When we do not consider neutrino oscillations, the yield is 1.50 x 1O~ 7 M . The yield of 7 Li 
increases by a factor of 1.19 owing to the neutrino oscillations. 

The mass fraction of 7 Be is larger by about a factor of 2.5 than that without the neutrino 
oscillations in the innermost region of the He layer, 3.8M < M r < 4.3M (see Fig. 3a). In 
the range of 4.3M < M r < 4.6M , where the 7 Be mass fraction increases with the mass 
coordinate, the increase in the 7 Be mass fraction due to the neutrino oscillations is more 
than a factor of 3. In the range of M r > 4.6M it is a factor of 2.5 ~ 3.2. Finally, the 
obtained yield of 7 Be is 2.66 x 1O~ 7 M . Without the neutrino oscillations, the yield of 7 Be is 
8.62 x 1O _8 M . The yield of 7 Be increases by a factor of 3.1 owing to the neutrino oscillations. 
Thus, the total yield of 7 Li, the sum of the yields of 7 Li and 7 Be, is 4.45 x 1O~ 7 M with 
the neutrino oscillations and 2.36 x 1O _7 M without the oscillations. The 7 Li yield with the 
neutrino oscillations is larger by a factor of 1.89 than that without the oscillations. 

The mass fraction of 11 B in the range of M r > 4.3M of the He layer is larger by about 
a factor of 1.2 ~ 1.3 than that without the oscillations (see Fig. 36). Inside the range of the 
He layer, it increases by about a factor of 2.1. The increase in the mass fraction of n C in the 
He layer is about a factor of 3.4 ~ 3.9. This increase is larger than that of n B and is close 
to the factor of 7 Be. The n B and n C mass fractions in the O-rich layers of M r < 3.8M are 
not affected strongly by the neutrino oscillations. The yields of n B and n C in this case are 
7.02 x 1O _7 M and 9.16 x 1O _8 M . Without the neutrino oscillations, the corresponding 
yields are 5.63 x 1O" 7 M and 6.29 x 1O _8 M . Thus, the total n B yield increases by a factor 
of 1.27. 

We discuss the production process of 7 Li and n B in the He layer and the effect of the 
neutrino oscillations on the //-process reactions. In the He layer, 7 Li and 7 Be are produced 
through 4 He(z/, u'p) 3 B.(a, 7) 7 Li and 4 He(z/, u'n) 3 B.e(a, 7) 7 Be, respectively. The corresponding 
charged- current z/-process reactions are 4 He(z/ e , e + n) 3 H(o;, 7) 7 Li and 4 He(z/ e , e~p) 3 He(o;, 7) 7 Be. 
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Most of n B is produced through 7 Li(o;,7) 11 B and the contribution from 12 C(z/, i/'p) n B is 
small. The n B in the O-rich layers is also produced from 12 C. The n C is produced through 
12 C(z/, u'n) u C The charged-current //-process reactions from 12 C are 12 C(z/ e , e + n) 11 B and 
12 C(z/ e , e~p) 11 C. Main production processes of these light elements are also written in Yoshida 
et al. (2004, 2005, 2006). In the case of normal mass hierarchy and adiabatic H-resonance for 
neutrinos, v e and v^ T are completely exchanged owing to the adiabatic mixing, and therefore 
the rates of 4 He(z/ e , e~p) 3 He and 12 C(z/ e , e~p) u C increase. Thus, the mass fractions of 7 Be 
and n C become larger. On the other hand, there are no resonances for antineutrinos. The 
transition probability between V e and v^ T is small so that the rates of 4 He(z/ e , e + n) 3 H and 
12 C(z/ e , e'p) u C scarcely become larger, and therefore the increase in the mass fractions of 
7 Li and n B is small. The increase by a factor of 2 in the mass fraction of n B is due to the 
production by way of 7 Be(n, p) 7 Li(a, 7) n B. 

Secondly, we consider the case of normal mass hierarchy and sin 2 29 1 3 = 1 x 10~ 6 (nona- 
diabatic H-resonance for neutrinos) shown by thin lines in Figs. 3a and 36. The increment 
degree of 7 Li mass fraction compared with the case without the neutrino oscillations be- 
comes gradually larger with increasing the mass coordinate (see Fig. 3a). The maximum 
increment degree is 1.3 at the outer edge of the He layer. The mass fraction of 7 Be is larger 
by a factor of 1.2 compared with the case without the neutrino oscillations where the mass 
fraction has the maximum value. The increment degree becomes larger with increasing the 
mass coordinate, but the mass fraction is much smaller there. Thus, the obtained yields of 
7 Li and 7 Be are 1.65 x 1O~ 7 M and 1.02 x 10~ 7 M Q . The total yield of 7 Li is 2.67 x 1O~ 7 M ; 
it is larger by a factor of only 1.13 than the case without the neutrino oscillations. 

The increment degree of the n B mass fraction due to the neutrino oscillations gradually 
increases with the mass coordinate in the He layer (see Fig. 36). However, increment degree 
is 1.2 at the maximum at the outer edge of the He layer. This is due to smaller transition 
probability from u^ T to u e (see Fig. 16). The dependence of increasing n C mass fraction on 
the mass coordinate is similar to that of n B. The maximum increment degree is 2.9. The 
obtained yields of n B and n C are 5.75 x 1O _7 M and 6.47 x 1O _8 M and, therefore, the 
total n B yield is 6.40 x 1O _7 M . It is only slightly larger by a factor of 1.02 compared with 
the case without the oscillations. 

Thirdly, we consider the case of inverted mass hierarchy and sin 2 29\% = 1 x 10~ 2 (adi- 
abatic H-resonance for antineutrinos) shown by thick lines in Figs. 3c and 3d. The mass 
fraction of 7 Li is larger by a factor of 1.4 ~ 1.5 in most region of the He layer. For n B, the 
mass fraction increases by a factor of 1.1 ~ 1.4. The yields of 7 Li and n B are 2.22 x 1O~ 7 M 
and 7.05 x 1O~ 7 M . Compared with those without the oscillations, these yields are larger 
by factors of 1.48 and 1.25. Although adiabatic H-resonance appears for antineutrinos and, 
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therefore, the flavor change between v e and v^ T occurs completely, the increment degree of 
the 7 Li and 11 B mass fractions are smaller than those of 7 Be and n C in the corresponding 
case of normal mass hierarchy. This is because the difference between the mean energies for 
v e and is smaller than that for v e and z/^ jT . This small difference of the mean neutrino 
energies causes smaller enhancement of 7 Li and n B even in the adiabatic mixing. 

The mass fractions of 7 Be and n C are larger than those without the neutrino oscillations 
in the range of M r > 4.8M in the He layer. On the other hand, they are smaller than 
those without the neutrino oscillations inside the region. This is because neutrons produced 
through 4 He(z/ e , e + n) 3 H, of which rate is enhanced by the neutrino oscillations, decompose 
3 He, 7 Be, and n C. The yields of 7 Be and n C are 8.20 x l(r 8 M and 6.35 x 1O" 8 M . The 
yield ratios to the case without the neutrino oscillations are 0.95 and 1.01 for 7 Be and n C. 
The total yields of 7 Li and n B are 3.04 x 1O~ 7 M and 7.69 x 1O~ 7 M , and their increment 
factors are 1.29 and 1.23. 

Finally, we show the case of inverted mass hierarchy and sin 2 2#i 3 = 1 x 10 -6 (nonadia- 
batic H-resonance for antineutrinos) represented by thin lines in Figs. 3c and 3d. The mass 
fraction distributions of all four nuclear species are the same as the corresponding distribu- 
tions in the case of normal mass hierarchy and nonadiabatic H-resonance. We showed in §3 
that the exchange probabilities for u e (u e ) — > f e (P e ) and u e (u e ) — > ^, T (^,r) are identical in 
the cases of nonadiabatic H-resonance independent of mass hierarchy at e v = 50 MeV. The 
same conclusion is inferred for neutrinos at different neutrino energies. 



4.2. Yields of 7 Li and n B 

In this subsection we discuss detailed dependence of the 7 Li and n B yields on the mass 
hierarchy and sin 2 29 13. When we do not consider neutrino oscillations, the calculated yields 
of 7 Li and n B are 2.36 x 10~ 7 and 6.26 x 1O~ 7 M . We evaluate the ratios of the 7 Li and n B 
yields with and without neutrino oscillations. 

Figure 4 shows the dependence of the yield ratios of 7 Li (a) and n B (b) on sin 2 26>i 3 in 
the range between 1 x 10~ 6 and 1 x 10 _1 . We find three characteristics of the 7 Li and n B 
yields as to the dependence on sin 2 2# 13 . In the case of sin 2 26* 13 < 2 x 10~ 5 , the 7 Li and 
11 B yield ratios keep constant values and do not depend on mass hierarchy. In the range of 
2 x 10~ 5 < sin 2 2# 13 < 2 x 10~ 3 , the 7 Li and n B yield ratios increase with sin 2 2# 13 . The 
difference due to mass hierarchy is also seen. In the range of sin 2 29± 3 >2x 10~ 3 , the yield 
ratios roughly keep constant values again. The difference due to mass hierarchy is also seen 
in this range. The above dependence is compared with that of the adiabaticity parameter 
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7 (see Eq. (10)) and the flip probability Pf (see eq. (11)) of H-resonance at the location 
of the resonance density. The value 1 — Pf related to sin 2 29 13 is also drawn in Fig. 4. We 
set the neutrino energy to be 50 MeV. This neutrino energy is close to the optimum energy 
contributing most strongly to the charged-current //-process reactions of 4 He. 

In the case of sin 2 2#i 3 <2x 10™ 5 , 1 — Pf is almost equal to 0; H-resonance is nonadi- 
abatic. In the case of 2 x 10~ 5 < sin 2 2# 13 <2x 10~ 3 , 1 — Pf increases with sin 2 20 13 , i.e., 
the resonance changes from nonadiabatic to adiabatic with increasing sin 2 20 13 . We call this 
range of sin 2 2# 13 "transition range" . In the case of sin 2 2# 13 >2x 10~ 3 , the value of 1 — Pf 
is almost equal to 1, i.e., H-resonance is adiabatic. The change of the flip probability Pf as 
a function of sin 2 2#i 3 is roughly similar to the change of the yields of 7 Li and n B. Thus, 
we conclude that the dependence of the 7 Li and 11 B yields on sin 2 2#i 3 strongly correlates to 
the adiabaticity of H-resonance. 

In the case of the nonadiabatic range, sin 2 26>i 3 <2x 10~ 5 , the 7 Li and n B yield ratios 
are about 1.1 and 1.02, respectively, and they are independent of mass hierarchy. In the 
limit of sin 2 20 13 = 0, the yield ratios still stay near the above values. We showed in §3 that 
the transition probabilities do not depend on mass hierarchy in nonadiabatic region. This 
characteristics can also be seen in the 7 Li and n B yield ratios. The yield ratios 1.1 and 1.02 
are slightly larger than unity, which reflect small enhancement of the mass fractions of 7 Li 
and n B in outer region of the He/C layer as shown in §4.1. 

In the transition range, the yield ratios of 7 Li and n B increase with sin 2 20 13 and a 
difference relating to mass hierarchy appears. The increase in the 7 Li yield is due to the 
enhancement of the 7 Be production through 4 He(z/ e , e~p) 3 He(a, 7) 7 Be. The increase in the 
n B yield also arises from the enhancement of 7 Li yield by way of 7 Be. 

In the adiabatic range, the 7 Li yield ratio depends on mass hierarchy: the 7 Li yield ratio 
is 1.9 in the case of normal mass hierarchy, and 1.3 in inverted mass hierarchy. As shown in 
§4.1, the increase in the reaction rate of 4 H(z/ e , e~p) 3 He raises the 7 Be production in normal 
mass hierarchy. The increase in the rate of 4 He(z/ e , e + n) 3 H also raises the 7 Li production 
in inverted mass hierarchy. The increase in the reaction rate of 4 He(z/ e , e~p) 3 He in normal 
mass hierarchy is larger than that of 4 He(z/ e , e + n) 3 H in inverted mass hierarchy. In the case 
of T Ve = 6 MeV, the reaction rate of 4 He(z/ e , e~p) 3 He is larger than that of 4 He(z/ e , e + n) 3 H. 
Further, the difference between T Ve and T V ^ T is larger than that between T Pe and T Pfi T (where 
note that T v = T p ). Thus, the increase in the 7 Li yield ratio in normal mass hierarchy 
is larger than that in inverted mass hierarchy. The n B yield ratio also depends on mass 
hierarchy, but the difference is smaller than that of the 7 Li yield ratio. As mentioned above, 
the increase in the 11 B yield arises in this case from enhanced 7 Li production by way of 7 Be. 
However, most of 7 Li produced through 7 Be do not capture a-particles to form n B. 
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4.3. Dependence on Initial v e and v e Temperatures 

In the present study we adopted the temperatures of v e , 9 e , and v^ T (P^ T ) equal to 
3.2 MeV, 5.0 MeV, and 6.0 MeV, respectively. On the other hand, some other studies 
(e.g. Woosley & Weaver 1995; Rauscher et al. 2002) adopted the temperature of u e and 
P e as T Ue = T Pe = 4.0 MeV. Since the enhancement of the 7 Li and n B yields depends 
on the neutrino temperatures at the neutrino sphere, we can find some effects on different 
temperatures of u e and v e . We consider four sets of the temperatures of v e and v e : (T Ue ,T Pe ) 
= (3.2 MeV, 5 MeV), (4 MeV, 4 MeV), (4 MeV, 5 MeV), (3.2 MeV, 4 MeV) as mentioned 
in §2.1. The first set is our standard model. The second set is the one adopted in Woosley 
& Weaver (1995), Rauscher et al. (2002), and so on. The third set satisfy T Ve < T Pe , but 
the difference of the temperature is smaller than the first set. The fourth set is prepared 
for comparison to the other three models: either T Ve or T Pe is different from the other sets. 
In these four sets we use a temperature of v^ T and v^ T , 6.0 MeV. The 7 Li and n B yields 
without neutrino oscillations in the four cases of the neutrino temperatures are listed in 
Table 2. 

Figure 5 shows the 7 Li and n B yield ratios as functions of sin 2 2#i 3 in the above four 
cases of neutrino temperatures. In the case of normal mass hierarchy (see Fig. 5a), the 7 Li 
yield ratio depends on initial temperatures in the nonadiabatic range, and the dependence 
becomes larger in the adiabatic range. In the case of small temperatures of v e and P e , i.e., 
(T„ e , T P J = (3.2 MeV, 4 MeV), the 7 Li yield ratio becomes 2.1 at the maximum. On the 
other hand, large neutrino temperature, i.e., (T Ue ,T Pe ) = (4 MeV, 5 MeV), provides the 7 Li 
yield ratio of 1.8. In the adiabatic range, the final 7 Li yield does not depend on the initial 
neutrino temperatures. The difference is due to the difference of the 7 Li yields without 
neutrino oscillations. After the neutrino mixing at H-resonance, the v e energy spectrum 
becomes FD distribution with T Ve = 6 MeV, independent of the initial T Ve . In addition, the 
yield produced originally as 7 Be is much larger than the 7 Li yield produced without neutrino 
oscillations. Therefore, the 7 Be yield does not depend strongly on the neutrino temperatures 
among the four sets. 

In the case of inverted mass hierarchy (see Fig. 56), the dependence of the 7 Li yield ratio 
on the initial neutrino temperatures is similar to that in normal mass hierarchy. The yield 
ratio is about 1.4 in the case of (T„ e , T Pe )=(3.2 MeV, 4 MeV) and in the adiabatic range. 
It is 1.26 in the case of (T Ue ,T Pe ) = (A MeV, 5 MeV). In the adiabatic range the produced 
amount of 7 Li through 4 He(z/ e , e~n) 3 H does not depend on the initial neutrino temperatures. 
The different neutrino temperatures lead to different 7 Li yields through 4 He(z/, u'p) 3 B., which 
makes virtually larger effect than the neutrino oscillations. 

Figure 5 c shows the cases of the 11 B yield ratios in the case of normal mass hierarchy. 
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The n B yield ratios are 1.31 and 1.23 in the cases of (T Ue , T P J=(3.2 MeV, 4 MeV) and (4 
MeV, 5 MeV), respectively, in adiabatic range. Since most of n B is produced by way of 7 Li, 
the dependence on the neutrino temperatures is similar to 7 Li. Some 11 B are produced by 
way of 7 Be through 7 Be(n,p) 7 Li(o;, 7) n B. Thus, the difference of n B yield due to different 
neutrino temperatures is smaller in the adiabatic case than in the nonadiabatic case. From 
the viewpoint of the yield ratios, the n B yield ratio is larger in the adiabatic range. 

In the case of inverted mass hierarchy (see Fig. 5(f), the n B yield ratios are 1.39 and 
1.21 in the cases of (T Ve , T P J=(3.2 MeV, 4 MeV) and (4 MeV, 5 MeV) in the adiabatic range. 
In the case of (T„ e , T Pe ) = (3.2 MeV, 5 MeV) the n B yield ratio is 1.35 in the adiabatic range. 
Thus, the n B yield ratios in the cases of T Ve = 3.2 MeV are larger than the maximum n B 
yield ratio in normal mass hierarchy. In these cases the 7 Li yield ratios are correspondingly 
similar to each other. This reflects the fact that 11 B is mainly produced by way of 3 H and 



When we calculated neutrino oscillations, we adopted the density profile of a presuper- 
nova and did not consider the shock propagation during supernova explosion. One of the 
reasons is that our hydrodynamical model assumes the inner boundary of the supernova 
ejecta as a mass cut and the density structure inside the mass cut is not considered. Indeed, 
numerical simulations of core collapse, core bounce, and explosion of surrounding materials 
have not been successful yet (Buras et al. 2003; Sumiyoshi et al. 2005). On the other hand, 
numerical simulations of the shock propagation from a neutrino sphere to the envelope have 
been studied using time-dependent inner boundary, i.e., the surface of the neutron star (e.g., 
Tomas et al. 2004). In this section we discuss the shock propagation effect on neutrino 
oscillations using the density profile of the supernova ejecta calculated by hydrodynamical 
model, and assuming a simple analytical density profile inside the mass cut. 

In order to discuss the shock propagation effect on the neutrino oscillations, we take 
simple analytical density structure inside the mass cut. First, we assume that the density 
drops as p oc r -3 from a neutrino sphere until it drops to the value at the mass cut. Then, 
the density stays a constant at the value for the location of the mass cut. The density profile 
in the present discussion is thus 



7 Li. 



5. 



Discussion 



5.1. Neutrino Oscillations in Shock Propagating Medium 




for r < r c 



(13) 
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and 

p(r) = p mcut for r c < r < r mcut , (14) 

where p mcu t and r mcut are the density and radius at the mass cut and r c is the radius where 
the density first drops to p = p mcut - As the shock wave moves outward, the density behind 
the shock front decreases and the density at the mass cut also drops. Thus, r c becomes 
larger as time passes by. 

We now study the influence of the neutrino oscillations with the shock propagation to 
nucleosynthesis of the light elements. Figure 6 shows the dependence of the 7 Li and 11 B 
yield ratios on sin 2 26*13. First, we consider the dependence in the normal mass hierarchy. 
We find in the case of sin 2 29 13 > lx 10~ 4 that the yield ratios of both 7 Li and n B are slightly 
smaller than the corresponding yield ratios in the presupernova density profile without shock 
propagation effect. The decrease in the yield ratios is about at most 4% in the case of 
sin 2 20i3 ~ 1 x 10~ 3 . The shock propagation effect of the neutrino oscillations is seen in the 
case of H-resonance close to adiabatic. 

Next, we consider the case of inverted mass hierarchy. Compared with the case of normal 
mass hierarchy, the decrease in the yield ratios are smaller; they are 1% level in the case of 
sin 2 20i3 ~ 1 x 10~ 3 . Therefore, we can conclude that in the inverted mass hierarchy not 
only the effect of neutrino oscillations itself but also the shock propagation effect is smaller 
than the case of normal mass hierarchy. 

The shock propagation effect on the neutrino oscillations is seen in 5 ~ 10 s when the 
shock front reaches the O/C layer or the inner region of the He/C layer. Meantime, the 
shock wave passes through the H-resonance region. When the shock wave passes through 
this region, the density becomes higher and the resonance region moves outward in the mass 
coordinate. At the same time, the density change at the shock front makes the adiabaticity of 
the resonance slightly smaller. After the shock front has passed the O/C layer, the location 
of the resonance in the mass coordinate goes inward owing to the density decrease by the 
expansion. The shock propagation effect is mainly seen in the inner region of the He layer, 
3.8M Q < M r < 4.3M Q . We have shown in the last section that 7 Be is produced in the region 
of M r < 4.8M and that 7 Li is produced outside this region. In the normal mass hierarchy 
the increase in the yield is mainly due to 7 Be production. Therefore, the shock propagation 
effect is seen more clearly in the normal mass hierarchy. 
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5.2. Dependence on Stellar Models 

We investigated the ^-process nucleosynthesis with neutrino oscillations using 14E1 
model, which is a progenitor model of SN 1987A in Shigeyama & Nomoto (1990). This 
stellar model corresponds to about 20 M at the zero-age main sequence (ZAMS) and has 
a He core of 6.0 M . Stars with the mass larger than ~ 12M are considered to evolve 
to form Fe core and become core-collapse supernovae at the end. Their internal structure 
indicates "onion" shell structure and their abundance distribution depends on their stellar 
mass at the ZAMS. On the other hand, the treatment of convection also affects the structure 
of presuper novae. We here discuss the influence of the internal structure to the //-process 
nucleosynthesis with neutrino oscillations. 

Detailed studies of massive star evolution (e.g., Nomoto & Hashimoto 1988; Hashimoto 
1995) indicated that the region of the O-rich layers in the presupernova stage increases in 
larger ZAMS stellar masses. Hashimoto (1995) indicated that 4, 8, and 16 M He star 
models have the O-rich layers of 0.7, 4.4, and 11.4 M and that the region of the He-rich 
layers has commonly about masses of ~ 2M . Stellar models including semi convection and 
convective overshooting mixing also have large O-rich layers (Woosley & Weaver 1988). If 
the region of the O-rich layers is much larger than that of the He/C layer, most of n B would 
be produced through the //-process from 12 C in the O-rich layers rather than from 4 He in 
the He/C layer. On the other hand, the density in the O/C layer scarcely depends on the 
stellar mass (e.g., Nomoto & Hashimoto 1988; Hashimoto 1995) and the density is close to 
the density of H-resonance. We showed that the effect of neutrino oscillations is not seen 
in the O-rich layers. Thus, we expect that the increase in the 11 B yield due to neutrino 
oscillations is very small if the region of the O-rich layers is larger than the He/C layer and 
n B is mainly produced from 12 C. 

The density profile of the He/C layer also may affect the increase in the 7 Li yield by 
neutrino oscillations. This is because 7 Li is produced as 7 Be in the inner region of the He/C 
layer and as 7 Li outside the region. We expect that stellar models which produce more 7 Be 
rather than 7 Li bring about larger increase in the 7 Li yield due to neutrino oscillations in 
adiabatic H-resonance and normal mass hierarchy. In our model, the yield of 7 Be increases 
by a factor of 3.1 at the maximum, whereas 7 Li produced through 3 H(o;,7) 7 Li increases by 
a factor of 1.2 in the normal mass hierarchy. When we consider adiabatic H-resonance and 
normal mass hierarchy, stellar models that produce 7 Li as 7 Be rather than 7 Li lead to the 7 Li 
yield more than twice as that without oscillations. On the other hand, if the contribution 
from 7 Li through 3 H(o;,7) 7 Li is large, the increase in the 7 Li yield due to the oscillation 
would be smaller in the normal mass hierarchy, and it would be larger in the inverted mass 
hierarchy. In the latter case, the 7 Li yield produced through 3 H(o;, 7) 7 Li increases by a factor 
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of 1.5 at the maximum. 

The shock propagation effect on the neutrino oscillations would depend on the stellar 
mass. Less massive stars have smaller O-rich region, so that it takes shorter time until the 
shock wave arrive at the H-resonance density region. After the shock arrival, the transition 
probabilities to the other flavors become small and the increase in the rates of charged- 
current //-process reactions also reduces. Thus, the increase in the 7 Li and n B yields is 
expected to be small for supernovae evolved from less massive stars. 

5.3. CP Phase Effect 

In this study, we assume that the CP phase 5 is equal to zero because the definite 
value of 5 has not been determined from neutrino experiments. Let us briefly discuss the 
influence of CP phase on the 7 Li and n B yields produced through the //-process. We consider 
the transition probabilities to v e or P e . Yokomakura et al. (2002) showed exact relations of 
neutrino transition probabilities in arbitrary matter profile. They showed that the transition 
probability of v e — > v e (p e — > P e ) does not depend on the CP phase 5. From their study, we 
also obtained that the sum of the transition probabilities of — > v e and v T — > v e (p^ — > v e 
and u T — > P e ) does not depend on the CP phase. Here, we note that and u T (P^ and P T ) 
emitted from a proto-neutron star have a same energy spectrum, i.e., the numbers of and 
v T (P^ and P T ) are same for a given neutrino energy. Since the numbers of and u T (P^, and 
P T ) for a given neutrino energy are same, the transition probability from or v T (P M or P r ) 
to v e (p e ) for a given neutrino energy does not depend on the CP phase, and therefore, the 
change of the energy spectrum of v e (p e ) does not depend on the CP phase, too. We do not 
expect any influences of the CP phase to the 7 Li and n B yields. 

5.4. Active-sterile neutrinos 

It is noted that the existence of the fourth sterile neutrinos may change the effect of 
supernova light element synthesis on neutrino oscillations. Fetter et al. (2003) considered 
additional sterile neutrinos with the mass squared difference of 1 ~ 100 eV 2 from electron 
neutrinos and investigated the effect of active-sterile neutrino conversion on r-process in 
neutrino-driven wind model. The Los Alamos liquid scintillator neutrino detector (LSND) 
experiment suggested that neutrino oscillations occur with mass squared difference in the 
range of Am 2 = 0.2 ~ 10 eV 2 (Aguilar et al. 2001). They argued that it is difficult to 
generate neutrons enough to activate the r-process. Then, they showed that active-sterile 
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neutrino conversion makes the winds favorable for the r-process even in the initial condition 
unfavorable to the r-process. This is because the active-sterile neutrino conversion ceases v e 
and, thus, the neutron depletion by n(u e ,e~)p is suppressed. 

If the active-sterile neutrino conversion occurs effectively just after the onset of super- 
nova explosion, then this conversion may reduce the number of v e and, therefore, it may 
reduce 7 Li and n B yields. It is expected that the active-sterile neutrino conversion occurs 
in inner region of the O-rich layer or even deeper region (p ~ 6 x 10 5 g cm~ 3 for Am 2 = 1 
eV 2 and e v = 20 MeV; see Eq. (10)). The decrease in the number of v e may reduce not only 
the yields of 7 Li and n B but also those of heavy neutron-deficient nuclei like 138 La and 180 Ta 
more effectively, which are produced through (z/ e , e~) reactions. In practice, the effect of the 
active-sterile neutrino conversion is still complicated because mixing angles of active-sterile 
neutrinos have not been determined and the conversion may strongly depend on the mixing 
angles. 



5.5. Prospects for Constraining Mass Hierarchy and # 13 

We obtained that the 7 Li yield is larger by about a factor of two by taking account 
of neutrino oscillations of normal mass hierarchy and adiabatic H-resonance compared with 
that without neutrino oscillations. If we can detect the 7 Li abundance or the abundance ratio 
of 7 Li to an element that is not affected by the //-process in stars which show clear traces of 
supernovae and if one compares them with the evaluated abundance or abundance ratio in 
supernova ejecta as we discussed, we may be able to constrain mass hierarchy and the mixing 
angle #13. We show the dependence of the number ratio of 7 Li/ 11 B on sin 2 2# 13 in Fig. 7. 
The 7 Li/ n B ratio should be better than the abundance of 7 Li itself for observations because 
the 7 Li/ 11 B ratio is rather insensitive to the uncertainties of supernova neutrinos (Yoshida et 
al. 2006). When we evaluated the number ratio of 7 Li/ 11 B, we considered the four sets of the 
temperatures of v e and u e discussed in §4.3. For all three ranges of normal and inverted mass 
hierarchies and without neutrino oscillations, the largest and smallest values correspond to 
the cases of (T„ e , T Ve ) = (3.2 MeV, 4.0 MeV) and (4.0 MeV, 5.0 MeV). The 7 Li/ n B ratio with 
(T Ue ,Tjy e ) = (4.0 MeV, 4.0 MeV) is larger than the result with the neutrino temperatures in 
our standard model. In the normal mass hierarchy and sin 2 2#i3 > 2 x 10~ 3 , the 7 Li/ 11 B ratio 
is larger than 0.83. The maximum values of the 7 Li/ n B ratio in the inverted mass hierarchy 
and without neutrino oscillations are 0.71 and 0.63. We note that the range of the 7 Li/ 11 B 
ratio deduced from the uncertainties of the v^ T temperatures, which have been investigated 
in Yoshida et al. (2006), is included in the range obtained in the present study. Therefore, the 
7 Li/ 11 B number ratio in the normal mass hierarchy and adiabatic H-resonance is larger than 
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that without the oscillations. The 7 Li/ n B ratio could be a tracer of normal mass hierarchy 
and relatively larger 9 13 , still satisfying the observational constraint of sin 2 2# 13 < 0.1, if the 
number ratio of 7 Li/ 11 B is precisely observed. 

Recently, measurements of B isotopic ratio in metal poor stars have been challenged by 
several groups (e.g., Rebull et al. 1998; Rebull, Primas, & Duncan 2000). Since supernovae 
provide large amount of n B with small amount of 10 B, the stars having traces of a supernova 
are expected to indicate 11 B/ 10 B ratio larger than that of the solar-system composition. 
Deficiency of B abundance with normal primordial 7 Li is reported in a metal poor star 
formed in an epoch when supernova nucleosynthesis dominated in the early Galaxy with 
quenched contribution from the cosmic ray interactions with interstellar medium (Primas, 
Duncan, & Thorburn 1998; Primas et al. 1999). The ratio of 7 Li/ n B in supernovae might 
be found also in primitive meteorites. 

Primitive meteorites should contain supernova-originating materials and 7 Li/ 11 B ratio 
might constrain the neutrino oscillation parameters. The solar-system 7 Li/ 11 B ratio is 3.11 
(Anders & Grevesse 1989). The 7 Li/ 11 B in the Galactic cosmic rays, which is one of the main 
production sites of Li and B, is evaluated to be 0.58 (Fields et al. 2000). These two values 
are out of our 7 Li/ 11 B range evaluated in the supernova //-process. Therefore, supernova- 
originating material should have 7 Li/ 11 B ratios different from those of the solar-system com- 
position and the Galactic cosmic rays, which is a favorable feature for our detectable 7 Li/ 11 B 
ratio in the range 0.6 < N( 7 Li)/N( n ~B) < 1.0 . Recently, presolar grains from supernovae 
have been found and isotopic ratios of C, N, O, Si, and Ti have been measured (e.g. Amari 
et al. 1992; Nittler et al. 1996). Measurement of B isotopic ratios in the grains has been 
attempted (Hoppe et al. 2001). Future isotopic and abundance measurements of Li and B 
would deduce 7 Li/ 11 B ratio in supernovae. 

We have to note, however, that there are many uncertainties to evaluate the 7 Li/ 11 B 
ratio and the Li yield in supernova ejecta. These values depend on the stellar mass at 
least theoretically. As discussed in §5.2, treatment of convection may also change stellar 
structure and abundance distribution. Aspherical explosion may produce different amount 
of light elements. In our study, the energy dependence of //-process reaction cross sections 
is very simplified: more precise evaluation of the energy dependence is required (Suzuki et 
al. 2006). In order to evaluate the enhancement factor due to neutrino oscillations, we need 
to construct more precise nucleosynthesis model in massive star evolution and supernova 
explosions. We also need to observe 7 Li/ n B ratio and Li abundance in stars having traces 
of a definite supernova. 

We showed the possibility of constraining parameters of neutrino oscillations from the 
viewpoint of nucleosynthesis. The enhancement of 7 Li/ 11 B ratio and 7 Li abundance would be 



- 21 - 



seen in the normal mass hierarchy and adiabatic H-resonance. There are different approaches 
to constrain neutrino oscillation parameters. One is to constrain neutrino masses from 
cosmological observations. Massive neutrinos produced in the early universe affect Cosmic 
Microwave Background (CMB) power spectrum and structure formation, i.e., the shape of 
matter power spectrum (review in Elgar0y & Lahav 2005). Constraining the total mass 
of mass eigenstates of neutrinos (m^ tot = mi + m 2 + m 3 ) from cosmological observations 
has been carried out. An upper limit of the total neutrino mass, m^ tot < 0.42 eV at 95% 
C.L., has been deduced by fitting cosmological parameters combined with the analyses of 
Wilkinson Microwave Anisotropy Probe (WMAP) CMB, galaxy clustering, galaxy bias, and 
Lya forest of the Sloan Digital Sky Survey (Seljak et al. 2005). In the near future higher 
precision cosmological observations would reduce upper limit of the total neutrino mass, 
from which the squared mass difference and mass hierarchy also would be determined. 

Observations of Diffuse Supernova Neutrino Background (DSNB) may constrain pa- 
rameters of neutrino oscillations, too. Although DSNB has not been detected yet, future 
improved neutrino detectors may detect DSNB. Yiiksel, Ando, & Beacom (2005) proposed 
a gadolinium-enhanced Super-Kamiokande detector by which average v e energy and the to- 
tal v e energy per supernova are to be measured after five years run. Further, Ando (2004) 
discussed that, in the inverted mass hierarchy and adiabatic H-resonance, the average v e 
energy is very different from the expected one without the neutrino oscillations. There- 
fore, combined with theoretical evaluation of the average v e and P M)T energies emitted from 
proto-neutron stars, the effect of the neutrino oscillations will be studied and the neutrino 
oscillation parameters will be constrained. 

We have these three different procedures to constrain neutrino oscillation parameters. 
First, the investigation of the ^-process nucleosynthesis in supernovae would provide a piece 
of evidence for normal mass hierarchy and adiabatic H-resonance. Second, DSNB mea- 
surement with improved neutrino detector would reveal a possible proof of inverted mass 
hierarchy and adiabatic H-resonance. Third, improved cosmological observations would clar- 
ify mass hierarchy. Thus, combining the three investigations with each other, we hope to 
constrain strictly the neutrino oscillation parameters in the future. 

6. Conclusions 

We studied light element nucleosynthesis through the //-process in supernovae taking 
account of neutrino oscillations. The parameters of neutrino oscillations were adopted from 
the evaluations through several neutrino experiments. We used a supernova model corre- 
sponding to SN 1987A and investigated the dependence of the 7 Li and n B yields on the 
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mixing angle 6>i 3 and mass hierarchy. The obtained results are summarized as follows: 

• Neutrino oscillations affect the yields of 7 Li and 11 B synthesized in supernova explo- 
sions. Especially, the 7 Li yield increases by a factor of 1.9 in the normal mass hierarchy 
and adiabatic H-resonance (sin 2 2# 13 >2x 10~ 3 ) compared with that without neutrino 
oscillations. The n B yield increases by a factor of 1.3. 

• In the inverted mass hierarchy, the increase in the 7 Li and n B yields is smaller: the 
yields of 7 Li and n B increase by factors of 1.3 and 1.2. 

• Neutrino oscillations in supernovae make the reaction rates of charged-current v- 
process reactions larger. The reaction rates of neutral-current z/-process reactions 
do not change. Thus, the final amounts of the //-process products increase by the 
neutrino oscillations. In our study, main important //-process reactions for the 7 Li and 
n B production are 4 He(z/, z/'j9) 3 H(o;, 7) 7 Li, 4 He(z/, i/n) 3 He(a, 7) 7 Be, 12 C(z/, u'p) u B, and 
12 C(V, u'n) u C. When we consider neutrino oscillations, the following charged-current 
//-process reactions become also important: 4 He(z/ e , e _ p) 3 He(o;, 7) 7 Li, 4 He(z/ e , e + n) 3 H 
(a, 7 ) 7 Be, 12 C(z/ e ,e-p) n C, 12 C(v e , e+n) n B. 

• The neutrino temperatures also affect the 7 Li and 11 B yields due to neutrino oscilla- 
tions. Large difference of the temperatures of v e and v^ T brings about larger increase 
in the yields compared with those without neutrino oscillations. 

• The shock propagation effect on the neutrino oscillations would slightly reduce the 
increment of the 7 Li and n B yields. In our model, most of neutrinos pass through the 
He/C layer before the shock wave arrives at the O/C layer, i.e., the resonance density 
region. When the shock wave is in the O-rich layers, the density change by the shock 
does not influence strongly neutrino oscillations. 
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Fig. 1. — Transition probabilities by neutrino oscillations in the normal mass hierarchy. The 
value of sin 2 2^13 is (a) 1 x 10 -2 and (b) 1 x 1CT 6 . Solid lines are the transition probabilities 
from v e to u e , P ee , and from u^ T to u e ,P XC . Dashed lines are those from v e to u e , P ee b, and 
from v^ T to 9 e , P xe b- Horizontal axis is the mass coordinates in units of M Q . The neutrino 
energy is set to be 50 MeV. 
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Fig. 3. — Mass fraction distribution of 7 Li and n B in the normal mass hierarchy ((a) and 
(b)) and inverted mass hierarchy ((c) and (d)) as a function of the mass coordinate in units 
of M Q . O-rich, O/C, He/C, and He/N indicate the O-rich, O/C, He/C, and He/N layers in 
the supernova ejecta. Thick lines and thin lines correspond to the mass fractions of 7 Li and 
its isobar 7 Be, respectively, in (a) and (c), and those of n B and its isobar 11 C, respectively, 
in (b) and (of). Solid lines and dashed lines correspond to the cases of sin 2 2^13 = 1 x 10~ 2 
and 1 x 10~ 6 . The mass fractions calculated without neutrino oscillations are indicated by 
dotted lines. 
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Fig. 3. — Continued. 
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Fig. 4. — Dependence of the yield ratios of (a) 7 Li and (b) n B on sin 2 2^ 13 . Solid line and 
dotted line indicate the yield ratios in the normal and inverted mass hierarchies. In panel 
(a), dashed line indicates the flip probability in the case of neutrino energy equal to 50 MeV 
as the formulation of 1 — Pf. 
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Fig. 5. — Yield ratios of 7 Li ((a) and (6)) and n B ((c) and (cf)) related to the temperatures of 
and v e emitted from a proton-neutron star. Panels (a) and (c) correspond to the normal 
mass hierarchy and panels (b) and (d) correspond to the inverted mass hierarchy. Solid, 
dashed, dash-dotted, and dotted lines indicate the yield ratios in the cases of (T Ue ,T p J = 
(3.2 MeV, 5 MeV), (4 MeV, 5 MeV), (3.2 MeV, 4 MeV), and (4 MeV, 4 Mev). 
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Fig. 5. — Continued. 
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Fig. 6. — Dependence of the (a) 7 Li and (b) n B yield ratios on sin 2 29 13 taking account of 
the shock propagation effect on neutrino oscillations. Solid lines and dotted lines correspond 
to the normal and inverted mass hierarchies. Thick lines are the same as in Fig. 4. Thin 
lines indicate the yield ratios taking into account the shock propagation effect. 
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Fig. 7. — Dependence of the 7 Li/ n B number ratio on sin 2 2#i3. The dark- and medium- 
shaded regions indicate the ranges of the 7 Li/ 11 B ratios in the normal and inverted mass hi- 
erarchies. The number ratio without the neutrino oscillations is indicated by the light-shaded 
region. These regions are obtained from the evaluations using four sets of the temperatures 
of v e and v e defined in §4.3. Solid, dotted, and dashed lines in each region correspond to the 
number ratio with our standard model of the neutrino temperatures (see §4.3). See text for 
more details. 
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Table 1: Coefficients and threshold energies of the cross sections for charged-current reactions 
on 4 He and 12 C. The evaluated cross sections are in units of 10~ 42 cm 2 . 



Reaction 


a 




b 


e th (MeV) 


4 He(z/ e ,e" 


"p) 3 Hc 


2.45027 x 10" 


-6 


4.03496 


19.795 


4 He(z/ e , e 


f n) 3 H 


2.37357 x 10" 


-5 


3.31394 


21.618 


12 C(z/ e ,e" 


p) u c 


3.73589 x 10" 


-6 


3.78171 


17.939 




n) n B 


2.76414 x 10" 


-5 


3.08738 


17.761 


12 C(z/ e ,e" 


7 ) 12 N 


2.45749 x 10" 


-3 


1.88899 


17.3381 


12 C(u e ,e" 


7) 12 B 


2.13890 x 10" 


-3 


1.72299 


14.3909 



Table 2: Yields of 7 Li and n B without neutrino oscillations for four sets of the adopted 
temperatures of v e and V e . The temperature of u^ T and P MjT is assumed commonly to be 6.0 
MeV (see §2.2). 



(T Ue , T Ve ) 


M( 7 Li)/M 


M( n B)/M 


(3.2 MeV, 5 MeV) 


2.36 x 10" 6 


6.26 x 10~ 6 


(4 MeV, 4 MeV) 


2.24 x 10~ 6 


5.77 x 10- 6 


(4 MeV, 5 MeV) 


2.50 x 10~ 6 


6.80 x 10" 6 


(3.2 MeV, 4 MeV) 


2.08 x 10" 6 


5.23 x lO" 6 



